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A three-qubit 13 C solid-state nuclear magnetic resonance (NMR) system for quantum informa- 
tion processing, based on the malonic acid molecule, is used to demonstrate high-fidelity universal 
quantum control via strongly modulating radio-frequency pulses. This control is achieved in the 
strong-coupling regime, in which the timescales of selective qubit addressing and of two-qubit in- 
teractions are comparable. State evolutions under the internal Hamiltonian in this regime are 
significantly more complex, in general, than those of typical liquid-state NMR systems. Moreover, 
the transformations generated by the strongly modulating pulses are shown to be robust against the 
types of ensemble inhomogeneity that dominate in the employed molecular crystal system. The sec- 
ondary focus of the paper is upon detailed characterization of the malonic acid system. The internal 
Hamiltonian of the qubits is determined through spectral simulation. A pseudopure state prepara- 
tion protocol is extended to make a precise measurement of the dephasing rate of a three-quantum 
coherence state under residual dipolar interactions. The spectrum of intermolecular 13 C- 13 C dipolar 
fields in the crystal is simulated, and the results compared with single-quantum dephasing data ob- 
tained using appropriate refocusing sequences. We conclude that solid-state NMR systems tailored 
for quantum information processing have excellent potential for extending the investigations begun 
in liquid-state systems to greater numbers of qubits. 

PACS numbers: 03.67.Lx,61.18.Fs,76.60.-k 



I. INTRODUCTION 



Quantum information processing (QIP) aims to 
achieve the ultimate in computational power from phys- 
ical systems by exploiting their quantum nature |l| . Nu- 
clear magnetic resonance (NMR)-based QIP has been 
successfully implemented in liquid-state ensemble sys- 
tems of up to 7 qubits @j 0,01111 El 0- Universal quantum 
control is achieved through application of external radio- 
frequency (RF) fields on or near resonance with spin tran- 
sitions of a set of separately addressable, coupled spins. 
State initialization (to a fiducial state such as |00. .0) ) 
is effectively achieved in these systems by the prepara- 
tion of pseudopure states 0, El ■ Pseudopure states have 
recently been demonstrated in a 12-qubit liquid system 
Q and in a 12-spin liquid-crystal system [tj. A hall- 
mark of control in liquid-state systems is a separation 
of timescales between (faster) qubit addressability and 
(slower) two-qubit coupling gates, typically by an order of 
magnitude or more For the homonuclear subsystems, 
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this corresponds to the weak-coupling regime, in which 
the relative Zeeman shifts in the qubit Hamiltonian are 
significantly larger than the J-couplings 0, 0] . In this 
regime, the evolution due to spin interactions is predom- 
inantly of the controlled-phase form. In this work, we 
examine NMR-based QIP implemented in a solid-state 
homonuclear system, in which the qubit Hamiltonian is 
no longer in the weak-coupling regime. A solid-state 
NMR architecture is attractive due to several key prop- 
erties 0: (1) nuclear spin states have been purified to 
near-unity polarizations in solids 13]; (2) intrinsic deco- 
herence times can be much longer, and two-qubit gate 
times much shorter, than those in the liquid state; (3) 
the qubit spins can be brought into well-controlled con- 
tact with a thermal-bath of bulk spins, enabling entropy- 
reducing operations such as algorithmic cooling 0, Gil 
and quantum error correction |l6t Il7j to be carried out. 
The system we will describe is specially tailored so that 
the ensemble description of the system is-to a good 
approximation-analogous to that of liquid-state NMR- 
QIP, and therefore the general aspects of control and 
measurement are the same 0. A similar three qubit 
system using single -cry stal glycine was first explored by 
Leskowitz et. al. 18], in which the homonuclear two- 
carbon system was approximately weakly-coupled. How- 
ever, we will demonstrate that universal, coherent control 
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can be implemented in the strong-coupling regime: the 
regime in which the timescales of qubit addressing and 
qubit coupling are comparable. In this regime, the trans- 
verse spin operator terms from qubit-qubit dipolar cou- 
plings are less suppressed by the relative Zeeman shifts. 
These residual 'flip-flop' terms crLcr_ + crLa*L (where 
<j± = <j x ± ia y and {a x ,a y } are Pauli matrices) render 
the state evolutions more complex, in general, than those 
of weakly-coupled s pin systems. We show that strongly 
modulating pulses [I3 succeed in controlling the solid- 
state QIP system, with generality and with high-fidelity. 
Moreover, the pulses provide significant robustness for 
the desired transformations against the ensemble inho- 
mogeneities that are typical of solid-state NMR systems. 
This latter property is of great importance to the prac- 
tical application of quantum algorithms in such systems, 
and we present here a first step in its study. This control 
methodology is a key ingredient in the realization of solid- 
state NMR-QIP testbed devices, but could also extend to 
many other potential systems for quantum information 
processing. We remark that liquid-crystalline NMR-QIP 
implementations 0, H3| represent a coupling regime that 
is typically intermediate between the strong- and weak- 
coupling cases. Strongly modulating pulses would there- 
fore also be an appropriate means for controlling such 
systems universally. 

The secondary focus of this paper is to characterize, 
in detail, the employed three-qubit system based on the 
malonic acid molecule. This includes characterization 
of the dominant ensemble inhomogeneities arising both 
through linear and bi-linear terms in the Hamiltonian. 
In tailoring the present system, we have used dilution of 
the qubit molecules as a means for reaching an approxi- 
mate ensemble description in which processor molecules 
arc, ideally, non-interacting and reside in identical envi- 
ronments. However, the need for a macroscopic number 
of spins to generate observable NMR signals through the 
usual inductive detection both limits the degree of di- 
lution and requires a macroscopic sample. The former 
results in perturbing intermolecular dipolar fields, and 
the latter typically yields significant distortions of the 
applied magnetic fields over the sample volume, namely, 
of the static magnetic field (due to susceptibility/shape 
effects) and of the RF amplitude (due to the RF coil 
geometry). In this paper, we address the robustness of 
strongly modulating pulse control to dispersion of Zee- 
man shifts and of RF amplitudes. We also characterize 
in detail the intermolecular dipolar environment in the 
present system, both theoretically and experimentally. 

The paper is organized as follows: section UTI reviews 
strongly modulating RF pulses as a means of achiev- 
ing universal quantum control, and numerical results 
for an example pulse are discussed; in section IIIII the 
dilute 13 C malonic acid system is first characterized; 
section IIVI demonstrates the preparation of a pseudop- 
ure state (as a benchmark for control) and analyzes 
the results; section IIV CI treats an application of the 
pseudopure state protocol: precise measurement of the 



triple-quantum dipolar dephasing rate, and comparison 
to single- quantum dephasing rates; in section [V] simu- 
lated and experimental data are presented that explore 
the effect of intermolecular dipolar couplings on qubit 
coherence times; in section IV Bl a multiple-pulse refo- 
cusing sequence is used in order to compare appropriate 
experimental quantities with the intermolecular dipolar 
simulations, and as a first step in testing the attainabil- 
ity of long ensemble coherence times in this system; the 
overall results are discussed in section IVT1 in the context 
of assessing the future goals and potential of solid-state 
NMR-QIP. 



II. UNIVERSAL CONTROL: STRONGLY 
MODULATING PULSES 

Numerically optimized 'strongly modulating' pulses 
were previously introduced as a means of implementing 
fast, high-fidelity unitary gate operations in the context 
of liquid-state NMR-QIP [l9|. The aim is to construct an 
arbitrary modulated RF waveform that, when applied to 
the system, generates a desired effective Hamiltonian cor- 
responding to a particular quantum gate. This is accom- 
plished numerically by minimizing the distance between 
the actual and the desired unitary transformations using 
a simplex search algorithm. Gate fidelities are defined by 
the expression: 

f = J2p^ u LuU)/n\ 2 (i) 

where N is the dimension of the Hilbert space, Udes is 
the desired unitary and U ca i c is the unitary calculated 
for the evolution of the system under the modulated RF 
pulse. Here, is a normalized empirical distribution 
over an inhomogeneous parameter (or parameters) of the 
ensemble, typically the RF amplitude. The expression 
for F corresponds to an average fidelity over all possible 
input states [l9| . The number of parameters the algo- 
rithm must search over is made minimal by requiring the 
modulating waveform to consist of a series of constant 
amplitude and frequency periods, so that each period 
has a time-independent Hamiltonian in a particular ro- 
tating reference frame . Modulation pulses with ideal 
(simulated) fidelities of order F ~ 99% are readily found 
for liquid-state NMR-QIP systems with up to 6 qubits 

Bill. 

It is observed empirically that good modulating pulse 
solutions tend to have the average RF amplitude vrf ~ 
\H\, where \H\ is the magnitude of the internal qubit 
Hamiltonian. It is precisely this strong driving regime 
in which analytical techniques (e.g. perturbation the- 
ory) for calculating dynamics break down, yet the system 
generally has the broadest (and most rapid) access to the 
manifold of allowed states. It is also the regime in which 
all accessible spin transitions are excited, so that refo- 
cusing of unwanted interactions becomes possible, even 
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to the extent that the dephasing effects of ensemble in- 
homogeneities may be partially suppressed [2(3, • 

An example of a strongly modulating RF pulse is 
shown in Fig. ^ The pulse is tailored to generate a three- 
qubit controlled- (not® not) gate in the strongly- coupled 
malonic acid system to be described in the next section. 
The average RF amplitude is 9.4 kHz, whereas the mag- 
nitude of the internal 13 C Hamiltonian (parameters listed 
in Table [I] next section) is \TL\ ~ 7.3 kHz. Part (c) of 
the figure demonstrates the robustness of the calculated 
unitary over the dominant inhomogeneous Hamiltonian 
parameters of the ensemble, namely, scaling of the RF 
amplitude and offset of the static field. This pulse was 
optimized over a 5-point probability distribution of 
RF amplitude scaling (corresponding to that measured 
in our RF coil), centered on unity, with standard devi- 
ation cr(p M ) ~ 6%. The fact that the unitary fidelity is 
> 90% over a 1 kHz range of static field offset demon- 
strates the ability of the modulated pulse to effectively 
rcfocus evolution under Zeeman shift dispersion (no Zee- 
man distribution was used in the optimization) . An ideal 
fidelity ~ 98% was obtained here, however, even greater 
precision is likely required for successful general imple- 
mentation of quantum algorithms. We consider this a 
first step that can be significantly improved upon in fu- 
ture work. 

The strongly modulating pulse methodology generates 
fast pulses, relative to traditional selective pulse meth- 
ods, to implement unitary gates. Although the example 
pulse of Fig. n has a duration ~ 700/zs, we have found 
pulse solutions for the same gate (fidelities > 90%) as 
short as ~ 450/zs. This compares well with the same 
gate carried out as two separate controlled-not gates im- 
plemented using standard low-amplitude selective ir/2 
pulses and 7r/2 controlled-phase evolutions; we estimate 
such a gate would require at least ~ 1.25 ms to imple- 
ment. Optimal control theory has been used by other 
workers to design numerical procedures for constructing 
near time optimal state transformation pulses p2^ . 
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FIG. 1: Example of a modulated RF pulse solution for gener- 
ating a three-qubit gate in the malonic acid system, (a) quan- 
tum network representation of the gate, which is a controlled- 
(not®not), with C m as the control bit and C\ and C2 as target 
bits, (b) the quadrature components of the RF modulation 
waveform, labelled 'real' and 'imaginary', (c) Simulated fi- 
delity of the transformation as a function of (left) scaling of 
the RF amplitude and (right) offset of the static magnetic 
field. 



III. CHARACTERIZATION I: ENSEMBLE 
HAMILTONIAN 

A. Dilute 13C-labeled single-crystal malonic acid 

The system under study is a single crystal of malonic 
acid grown from aqueous solution with a dilute fraction 
of fully 13 C-labelled molecules. The main crystal used 
herein has dimensions of 4x1.5x1.5 mm 3 and a labelled 
molecule fraction of 3.2% (a similar crystal with f .6% di- 
lution was used in the experiments of section IV Bp . For 
the malonic acid molecule in the solid, the three carbon 
nuclei have distinct anisotropic chemical shielding ten- 
sors so that when placed in a large static magnetic field, 
crystal orientations can be found for which each carbon 
is separately addressable in frequency. Protons, of which 
there are four per molecule (and are 100% abundant in 



the crystal), can be used to cross-polarize the carbon 
spins, and are otherwise decoupled from the 13 C system. 
Figure [21 shows the malonic acid unit cell as determined 
by x-ray diffraction ,24|. The space group is P-f so that 
the two molecules in the unit cell are related by inver- 
sion symmetry, and are therefore magnetically equiva- 
lent. The crystal orientation is chosen to maximize the 
intramolecular 13 C- 13 C dipolar couplings and the rela- 
tive 13 C Zeeman shifts. The full spin Hamiltonian of the 
system is 

Hit) = H c + n H + Hen + H RF (i) , (2) 

where He and Tin are the 13 C and 1 H Hamiltoni- 
ans, TicH is the interspecies coupling Hamiltonian, and 
Hftir(i) is the time-dependent Hamiltonian of the exter- 
nal radio-frequency field. In many experiments, Tiupif) 
includes a strong field resonant with the *H spins so that 
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FIG. 2: Unit cell of malonic acid, with crystal axes indicated. 
The two molecules are related by inversion symmetry, and 
are therefore magnetically equivalent (P-l space group). The 
light (dark) atoms are carbon (oxygen), and the hydrogen 
atoms are not shown. A = 5.156A, B — 5. 341 A and C = 
8.407101. 
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TLch is effectively removed from the Hamiltonian. The 
13 C Hamiltonian can be decomposed as 

He = ricz + ricD + "cd > \ 6 ) 

where the terms on the right side, from left to right, are 
the Zeeman, the intramolecular dipolar, and the inter- 
molecular dipolar terms. The Zeeman and intramolecular 
dipolar terms dominate the natural 13 C Hamiltonian in 
the dilute 13 C crystal, and will be used (along with the 
RF) in the construction of quantum gates. The inter- 
molecular couplings act as perturbations on these terms, 
and their effects will be examined in section V. Denoting 
single-spin Pauli matrices as X = a x ,Y — <y y , Z = 
a z , the Zeeman and intramolecular dipolar terms may be 
expressed as 

'X i CD a = — (2Z m Z n - Y m Y n - X m X n ), (5) 

m<n<3 

where Vj are the rotating-frame Zeeman frequencies and 
d mn are the dipolar coupling strengths. Table I lists these 
parameters, as obtained from fitting the 13 C spectrum of 
Figure |3 for the crystal orientation used throughout. It 
also lists the free-induction dephasing times, ; the cor- 
responding rates (T^*) -1 provide a measure of the degree 
of ensemble inhomogeneity. It will be seen in section V. 
that the dominant contribution to these rates is Zeeman 
shift dispersion. 



B. Experimental setup 

All experiments were carried out at room temperature 
on a Bruker Avance NMR spectrometer operating at a 



FIG. 3: 13 C-NMR spectrum (dotted line) and spectral fitting 
results (solid line) for 3.2% 13 C-labelled malonic acid in the 
crystal orientation used throughout. The fitting procedure 
adjusts the simulated Hamiltonian to match the spectrum, 
and uses Lorentzian line broadening to fit the (T^)^ 1 de- 
pashing rates. Proton decoupling was applied using a TPPM 
sequence |25^l at an RF power of 250 kHz. Note the most 
intense peaks, in the center of each multiplet, correspond to 
the natural abundance (1.1%) 13 C. The uneven heights of 
the labelled-molecule peaks indicate that the relative Zeeman 
shifts are insufficient to fully truncate the strong form of the 
intramolecular dipolar couplings; hence, the eigenstates of the 
system are near, but not equal to, the usual computational 
basis states. 

TABLE I: Zeeman shifts and intramolecular dipolar cou- 
plings of the 13 C Hamiltonian, columns 1-3 from left; dipolar 
couplings involving the methylene protons, columns 4-5; 13 C 
free- induct ion dephasing times (T|) and spin- lattice relax- 
ation times (Ti). Values in columns 1-5 are in kHz. 
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field of 7.1 T, and home-built dual-channel RF probe- 
head. The sample coil had an inner diameter of 3 mm, 
and the typical 7r/2 'hard' pulse lengths were 1.25^s and 
0.75^s for carbon and hydrogen, respectively. 

IV. PSEUDOPURE STATE 

A. Preparation method 

A labelled pseudopure state 0, El can be prepared 
in this system using a combination of tailored RF pulses 
and RF phase cycling (temporal averaging) . A schematic 
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FIG. 4: Schematic pulse sequence for pseudopure state prepa- 
ration using strongly modulating pulses. The states of the 
qubits are indicated in product operator format after each 
transformation. The state l 1 ! 2 ^*" is denoted as 'IIZ', for ex- 
ample. Strongly modulating pulses are represented by wave- 
forms labelled 'IIZ/ZZZ' and 'TQPP', whose durations are 
0.54 ms and 1.10 ms, respectively. A collective tt/2 rotation 
about the rotating-frame y-axis is labelled 'Y90'. The angle 
cj> represents a six-fold phase-shifting of the RF to filter out 
all coherences except the 3Q coherence. r(3Q) indicates a 
variable delay subsequent to the creation of the 3Q state that 
is used to measure that state's decay. 



of the pulse sequence is shown in Figure^] The first step 
is a selective transfer of 1 H polarization to the methylene 
carbon C m , utilizing the 1 H/ 13 C coupling Hamiltonian 
T~Lch- While it is not necessary for this protocol to begin 
with a selective transfer, we demonstrate it here because 
it may be used in other algorithms for controlled, selec- 
tive coupling of the qubit system to the bulk 1 H system 
(which can be considered as a thermal-bath of spin po- 
larization). To accomplish the transfer, a pulse sequence 
is applied on both 1 H and 13 C channels synchronously, 
which, by design, effectively removes Tic + Ti-H from 
the ensemble Hamiltonian |26l |28| . In such sequences, 
homonuclear dipolar terms are refocussed by toggling the 
interaction-frame Hamiltonian along the rotating-frame 
x, y and z directions, spending equal time along each 
axis. Therefore, the dual sequence creates (to lowest or- 
der in the Magnus expansion of the average Hamiltonian 
pjf) an effective 1 H- 13 C exchange Hamiltonian 



E 



d jk (Z^Z k + Y' ] Y k + XiX k 



(6) 



where djk are pairwise 1 H- 13 C dipolar coupling con- 
stants, and the indices j,k run over all 13 C, 1 H spins, 
respectively. In the special case that there is only one 
coupled carbon/proton pair with a coupling of dcu , ap- 
plication of the sequence for a time r = 3 /(Aden) will re- 
sult in a state exchange between the two nuclei. Since this 
is approximately the case for the strong methylene X H- 
13 C coupling in our oriented malonic acid system (see Ta- 
ble I), we may implement a nearly selective polarization 
transfer to C m of magnitude P ~ Pjj — 3.98Pc (Pc and 
Ph are the thermal equilibrium carbon and proton polar- 
izations, respectively). Furthermore, this selective trans- 



fer removes a very small amount of polarization from the 
*H bath, since only the methylene protons on a dilute 
fraction of 13 C-labelled molecules lose their polarization. 
The remaining bulk *H polarization is preserved since 
the sequence is effectively a time-suspension sequence for 
the bulk spins. In our system, a selective polarization 
transfer of about 83% of Pjj to C m was achieved with a 
12-p ulse subsequence [3(]] of the Cory 48-pulse sequence 
|26| : the duration of the 12-pulse sequence was 40 ps for 
maximum transfer. The thermal equilibrium 13 C polar- 
ization can be removed prior to such a transfer by rotat- 
ing the equilibrium 13 C state (oc Z) into the transverse 
(x — y) plane and allowing it to dephase under local *H 
dipolar fields. 

Product operator terms denoting qubit states are or- 
dered as C\ ® C% ® C m ; for example, the symbol \XZ 
corresponds to the state ll 1 ® X 2 ® Z m . Also, single-spin 
terms such as X 1 <g> ll 2 <g> l m are sometimes abbreviated 
as X , for example. The polarization transfer is followed 
by a 13 C modulating RF pulse that transforms the state 
HZ to ZZZ, and then a collective 7r/2 pulse that ro- 
tates this to XX X. In addition to single-quantum (1Q) 
terms, this state contains the triple-quantum (3Q) co- 
herence rr+ o"+er+ + <7_<7_(7_, where tr + = X + iY and 
er_ = X — iY. A subsequent modulating pulse (de- 
noted 'TQPP' for 'triple-quantum to pseudopure') trans- 
forms the 3Q coherence into the labelled pseudopure 
state X(i+Z)(i+Z). This state is observable as a single 
NMR transition of G\ [3l|]. In order to cancel all other 
signals arising from the 1Q terms in the state XXX, 
phase cycling is used which exploits the n-proportional 
phase acquisition of an n-quantum state under z-axis ro- 
tation. Choosing desired coherence pathways using phase 
cycling is widely practiced in modern NMR spectroscopy 
[llj. By shifting the phase of the RF by <j> during the 
TQPP pulse, the unitary transformation generated by 
the pulse is rotated by <j> about Z' = Z 1 + Z 2 + Z m : 



Utqpp{4>) = RfUtqppRz 



(7) 



where Rf = e~ l( - z ~>^/ 2 . Note this is only true since the 
internal Hamiltonian of the system commutes with Z' , 
hence a z-axis rotation can be accomplished by acting 
only on the phase of the RF pulse. We may decompose 
the state prior to the TQPP pulse into 3Q and 1Q terms, 
e.g. XXX = p±3Q + p±iq. The final density matrix is 
calculated as 

Pf{4>) = RtU tqP pR7' i '(p±3Q + p±iQ)Riu} qpp R~ <t ' 
= Rf ((e- 3i *<r+ + e 3i V_)(ll + Z)(J +Z)/8 + e^p')TCj, 

(8) 

where p' = U tqpp {p±iQ)Ul qpp . Incrementing <f> in units 
of | for 6 scans, alternately adding and subtracting, 
adds constructively the 3Q terms while cancelling the 
1Q terms, since 
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(-l)*- l e ±mfc 3 = 6 • S(n - 3m), 



(9) 
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FIG. 5: Experimental results for the pseudopure state pro- 
tocol, (a) readout of the state UZ just after the selective 
polarization transfer, with readout consisting of a swap from 
Cm to Ci and collective n/2 pulse about the y-axis to produce 
JOl; (b) readout of the state ZZZ by a Ci-selective tt/2 ro- 
tation to produce the observable state XZZ; (c) the labelled 
pseudopure state yields a single absorption peak from the C\ 
multiplet. Fitting the data by spectral simulation indicates a 
protocol state-correlation of 87% ± 5%, using spectrum (a) as 
a reference. 

where m is an odd integer, and here m = 1 since n < 3. 
The remaining z-rotation Rf is undone by incrementing 
the phase of the receiver along with that of the TQPP 
modulating RF pulse. The labelled pseudopure state is 
thus obtained as 

1 6 

-^(-l) fc - 1 i?- fc7r/3 p / (fc 7 r/3)i?f /3 = X(l+Z)(l+Z)/4:. 

k=l 

( 10 ) 

The amount of signal contained in the resulting pseu- 
dopure state is 1/4 that of the input state, equal to the 
proportion of the 3Q part of the state XXX. 

B. Pseudopure state results 

Results for the pseudopure state protocol of Fig. 0] 
are shown in Fig. [SJ It shows (a) readout of the state 
WZ just after the selective polarization transfer, with 
readout consisting of a state-swap from C m to C\ 
followed by a collective tt/2 pulse to produce 111; (b) 
readout of the state ZZZ by a C\ -selective 7r/2 rotation 
to produce the observable state XZZ] (c) the labelled 
pseudopure state yields a single absorption peak from 
the Ci multiplet. Using the spectrum (a) as reference, 
the state-correlation of the pseudopure state protocol, 
determined by spectral fitting, is 87% ± 5%. Similarly, 
the state-correlation of the TQPP pulse, using the 
XZZ spectrum (b) as a reference, is 97% ± 5%. These 




milliseconds 

FIG. 6: Triple-quantum (3Q) and single-quantum dephasing 
data, plotted as the natural log of the signal intensity. Expo- 
nential fits indicate decay time constants (in milliseconds) of 
[2.37 (3Q), 5.37 (C m ), 9.07 (C 2 ), 8.66(d)]. 

results serve as a benchmark for quantum control in the 
strongly-coupled dilute molecular crystal system. They 
suggest that dephasing due to Zeeman shift dispersion 
is largely suppressed by these pulses, since the average 
free-induction dephasing time of the qubits is T 2 * ~ 2 
ms, and the total duration of the two modulating pulses, 
~ 1.6 ms, is a large fraction of that time. 



C. Dipolar dephasing of triple-quantum coherence 

Dephasing of the 3Q coherence state a+a + a + + 
O-O-O- was measured by inserting a variable time 
delay r(3Q) following the 'Y90' pulse in Fig.0] A collec- 
tive 7r pulse was also inserted in the center of r(3Q) to 
refocus Zeeman Hamiltonian terms (Hahn echo), so that 
the measured dephasing rate reflects only the perturbing 
dipolar fields experienced by the nuclei. This 3Q state is 
significant from a QIP perspective since it consists of the 
two most fragile density matrix elements in the three- 
qubit system, the extreme off-diagonal elements that 
are contained in the so-called "cat-state" 1 000) + 1 1 1 1) . 
The resulting signal decay is shown in Fig. along with 
the 1Q dephasing data for each qubit, also measured via 
Hahn echo. The 1Q decay data were measured by first 
preparing the states X 1 , X 2 and X m , and using the 
pair of outer spectral peaks from each spin multiplet to 
gauge the signal, observing only at delays corresponding 
to multiples of the dipolar oscillation periods of these 
peaks. Exponential fits yield decay time constants (in 
milliseconds) [2.37 (3Q), 5.37 (C m ), 9.07 (C 2 ), 8.66(d)], 
with experimental uncertainty ~ ±11% for each value. 
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FIG. 7: Simulation of line-shapes resulting from intermolec- 
ular dipolar couplings in 1.9% molecular dilution 13 C mal- 
onic acid. The simulated data (solid lines) are shown with 
Lorentzian fits (dashed lines). The resulting linewidths 
(FWHM) for Ci, C 2 and C m are 9, 11 and 8 Hz respectively. 
Correcting for the larger spin concentration due to presence 
of natural abundance spins, the linewidths should be approx- 
imately 60% broader. 



Within error, the observed 3Q dephasing rate is equal to 
the sum of the three 1Q rates. In the next section, we 
will see that intermolecular 13 C- 13 C dipolar fields lead 
to an approximately Lorentzian frequency broadening 
of the 1Q transitions. In this regime, the 1Q dephasing 
appears as a Markovian process, and therefore we do 
not expect the 3Q rate to carry any information about 
correlations in the 3Q dephasing. The fact that the 
observed 3Q rate is the sum of the 1Q rates is consistent 
with this picture. 

Finally, we remark that the faster 1Q dephasing of the 
methylene carbon C m evident in Fig. is probably due 
to a residual interaction with its neighboring proton, as 
the methylene proton pair are strongly coupled (22 kHz) 
which makes it difficult to remove the 19 kHz C m -H m i 
coupling under the standard TPPM decoupling used 
here. This effect is directly evident in the data presented 
in section IV Bl (note, however, the data in section IV Bl 
was obtained in a different sample at a slightly different 
orientation with respect to the external field). 



V. CHARACTERIZATION II: 
INTERMOLECULAR DIPOLAR FIELDS 

A. Modelling intermolecular dipolar dephasing 

We now turn our attention to characterizing 13 C- 13 C 
intermolecular dipolar interactions in the dilute malonic 



acid system. The general Hamiltonian of the system is 

j<k i 

where djk — djk{rjk,Qjk) are the intermolecular dipolar 
coupling constants that depend on the internuclear vector 
of length rjk and orientation 6jk (with respect to the 
external field), hjk is the spin-operator of the form of 
the dipolar coupling as in Eq. and v% are the chemical 
shifts. Taking a reference spin I = 1 to be a C\ spin, 
we may transform to the C\ rotating frame by setting 
all vi = for I 6 {Ci}, and setting vy to the correct 
offset frequencies for V e {(72,63}, respectively. The 
latter offset frequencies vy are on the order of a few kHz 
(see Table [J), whereas the intermolecular couplings djk 
are much weaker (< 100 Hz) due to the rj k scaling of the 
dipolar interaction. We may therefore truncate couplings 
between unlike spins C a -C^ to the heteronuclear form, 
and rewrite the dipolar terms as an ensemble of reference- 
spin Hamiltonians (e.g. for the C\ reference ensemble): 

j^jk _ dajk yijya,t _|_ djjjk glj gPjk 

+ !hl!i(2Z lj Z' ak - X lj X^ jk - Y lj Y l: > k ), (12) 

where j indicates a particular reference spin from the C\ 
ensemble, and k denotes the set of spin-labelled molecules 
that interact with the j reference spin, a, (3 and 7 in- 
dicate spins belonging to the C m , C2 and C\ ensembles, 
respectively. Note that the above Hamiltonian does not 
include all dipolar terms in Eq. 1111 since it only includes 
strong-coupling terms between C\ spins, and so implicitly 
neglects much of the spin-diffusion dynamics. We wish to 
make an estimate of the dephasing rate (i.e. line broad- 
ening) of the reference spins due to this ensemble interac- 
tion Hamiltonian (see [23 for seminal work along these 
lines). An approximate result can be obtained in the 
spin-dilute regime by assuming that the reference spin 
only interacts with one nearby spin-labelled molecule, so 
that there is only one k value and = r H? div - 

Under this restricted model, we will now describe the 
dephasing process in terms of the quantum evolution of 
the ensemble. The system consists of the reference spin, 
denoted £, and the three interacting spins a,/3, 7. The 
initial state is described by the density matrix ® p a( a 7 , 
where at thermal equilibrium p a j3~ f = 11 is appropriate for 
high-temperature NMR. The unitary operator acting on 
the j th ensemble member of the system at time t is 

Uj(t) = e - m ^\ (13) 

The operators that act on the reduced density matrix 
of the reference spin, i.e. the Kraus operators, are 
derived from the submatrices of U: 

u mn = (^\u\r aM ), (14) 
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where iV'Sa-y) i s t ne n eigenvector of the ot,/3, 7 system where the number of interacting spins is q — 3 in our 



1/37 

in some eigenbasis. Since /0 Q ^ 7 = 1, the Kraus operators 
A rn are obtained by summing with equal weights over the 
basis vectors: 



n=l 



(15) 



case. In the standard computational basis, \tpai3j) 
{|000), |001), |111)}, we obtain for U 3 {t) operators of 
the form (subscripts fc, I, m S {0, 1} denote the basis vec- 
tor of the a, (3, 7 system): 
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|i)<o|<w 



L (-i)"+»%(-i)'-»Y( z£: 2 + e 2 

\ 1 



|0)(H^ ie ^(-^+^(-^-^( - e ' + e " ' )] (16) 



where = 2irtd n j, and <5 a b is the Kronecker delta. 

To study dephasing behavior, we apply the eight Kraus 
operators Ai, (t) to the reference spin state p% = X = 
|0)(1| + 1 1)(0|, obtaining 

4 =Y, A iimPd4iJ = £«jo><ii + %J1><0|), 

(17) 



kin I 



klm 



where 



,,'0 



iH*+»««(-i) . I_r 

2 

(18) 



and 



&i, =-x 

-«M-l)»-<fl w (-l)' e-^((-^-i/2) + e -^((-i)- + i/2) 

2 

(19) 

Inspection of equations El an d El shows that gains 
phases ±8 a j±9fSj±9jj/2 and ±9 a j±90j±39yj/2 charac- 
teristic of generic binomial distributions. Averaging over 
the ensemble of Hamiltonians T~t J dip -, we obtain p^(t) — 

Y^jLi pi{t) / N ■ Defining a correlation function F x (t) = 
Tr(p^(t) ■ X), its frequency spectrum is given by the 



Fourier transformation F x (lo) — f_ e wt F x (t). Equa- 
tions El an d El make clear that F x (ui) will simply reflect 
the ensemble distribution of the intermolecular coupling 
constants leading to frequencies 2n{±d a j ± dpj ± d-yj/2} 
and 2%{±d a j ± dpj ± 3d 73 /2}. 

To make a concrete calculation of the distribution of 
coupling constants d^j , let us further assume that the in- 
teracting molecule lies within the first shell of 26 neigh- 
boring unit cells. Each unit cell contains two (magneti- 
cally equivalent) molecules, giving a total of 156 atomic 



sites (159 after adding the 3 atomic sites of the unit cell 
neighbor to the reference molecule). We calculated the 
coupling constants to each of these sites using the Carte- 
sian atomic coordinates obtained by x-ray diffraction |24j 
and the unit cell vectors (A, B, C) (see Fig- EJ) - The cou- 
plings are explicitly of the form: 



dk{rk,9k) = 7c^ 



1 - 3cos 2 



(20) 



where 7c is the gyromagnetic ratio of 13 C, is the 
length of the intcrnuclear vector connecting spin k to 
the reference spin, and 8k is the angle between this 
vector and the external magnetic field direction. Note 
that since there are 53 molecular sites we are consid- 
ering, the random occupation of one site corresponds 
to a labelled-molecule concentration of = 1.9%, 
which is in the range of our sample concentrations. 
Following the discussion above, we constructed fre- 
quency histograms by averaging the spectral frequencies 
2ir{±d a j ± d/3j ± 3d 7 j/2} over the distribution of 
coupling constants to the 53 molecular sites. The 
same procedure was carried out for C%, Ci and C m as 
reference spins. Figure shows the resulting histograms 
calculated for the crystal orientation used throughout. 
Lorentzian functions provide good fits for the purpose 
of determining the spectral linewidths, as expected for a 
dilute spin system [12j. 

The simulated spectra of figure [7| correspond to the 
expected indirect-dimension line-shapes that would 
result in a two-dimensional NMR experiment under 
a Hahn echo sequence 11]. The fits yield full-width 
half-maximum (FWHM) linewidths Av ~ 10 Hz. The 
presence of sing le 13 C spins due to natural abundance 
(1.1% atomic concentration) adds to the total spin 
concentration. This broadens the linewidth estimate 
by a factor ~ 71 for a dilute labelled-molecule 
percentage rj < 10%, since linewidth is approximately 
linear in spin concentration in this dilute regime [12^ . To 
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B. Coherence times under dipolar+Hahn 
refocusing 




5 10 15 20 25 30 

milliseconds 



FIG. 8: Decay of signal intensity under the dipolar+Hahn 
refocusing sequence for natural abundance spins in the 1.6% 
isotopic dilution crystal. C m (a) and C m (6) had TPPM pro- 
ton decoupling RF power levels of 270 kHz and 330 kHz, 
respectively. Exponential fits to the data indicate effective 
coherence times (in milliseconds) of [150 (Ci), 120 (C2), 16 
(C m (o)), 30 (C m (b))], although clearly the values obtained for 
the slower Ci and Ci decays are more approximate than those 
for Cm- Dephasing times > 21 ms were not explored here due 
to the heating effects of high-power proton decoupling on the 
room temperature probe/sample. 



account for interactions with molecules beyond the first 
neighboring unit cells in this static broadening picture 
(i.e. ignoring spin-diffusion dynamics), we can make a 
crude shell-model approximation. The broadening from 
the n spins in a spherical shell of radius R, n(R) oc R 2 , 
will go as i? -3 \J n(R) oc R~ 2 . This is true since the 
width (i.e. standard deviation) of a binomial distribu- 
tion scales as y/n and the dipolar interaction scales as 
i?~ 3 . Hence, the linewidths should be larger by a factor 
YlkLi k~ 2 = 7r 2 /6. The results, adjusted for 7/ = 1.6%, 
are summarized in Table HT1 along with experimental data 
obtained in a 1.6% labelled-molecule dilution crystal. 
The experimental values indicate effective Lorcntzian 
linewidths of the natural abundance 13 C spins mea- 
sured via Hahn echo, and by a multiple-pulse sequence 
consisting of the MREV-8 dipolar refocusing sequence 
[2^| in combination with the Hahn echo (detailed in the 
following section). In the table, 'Simulation II' refers to 
the aforementioned spectral simulations. These values 
should correspond to the Hahn echo experimental data. 
'Simulation I' refers to simulations in which like-spin 
C a -C a couplings were omitted. These values should 
correspond roughly with the dipolar+Hahn refocusing 
data, as described in the next section. 



Application of a Zeeman and homonuclear dipolar re- 
focusing sequence was investigated, to compare with the 
simulations of the previous section, and to explore the 
attainability of longer 13 C coherence times. The em- 
ployed sequence is the MREV-8 [23 dipolar decoupling 
sequence in combination with a Hahn echo. The pulse 
spacing of the MREV-8 sequence was set such that one 
eight-pulse cycle was 240/^s in duration. A single 7r-pulse 
was applied in the center of a given dephasing period 
to refocus the effective Zeeman field [2^ of the MREV- 
8 sequence (the effective field is in the z — x plane, so 
that a 7r-pulse about the y-axis will refocus it). It is a 
quasi-evolution-suspension sequence because it refocuses 
like-spin dipolar evolution on a relatively short time-scale 
compared to the refocusing of Zeeman evolution. There- 
fore it does not fully remove effective dipolar couplings 
resulting from unlike-spin C a -Cp couplings (their mag- 
nitude will be substantially scaled, however). The data 
are shown in Figure [5] The decay of the magnetization 
of the natural abundance spins is shown, since the aim 
here is to see the effect of refocusing the m<ermolecular 
dipolar couplings on coherence times. Monitoring the 
natural abundance spins is ideal in this regard, since the 
average intermolecular dipolar environment is identical 
for all 13 C spins in the sample, and for the natural abun- 
dance spins, it is the only 13 C dipolar environment. On 
the other hand, efficient refocusing of the intramolecular 
couplings would require a shorter MREV-8 sequence, and 
therefore more pulses for a given dephasing time, yield- 
ing greater loss of signal due to pulse imperfections. The 
ability to suspend the evolution of the qubit spins (e.g. 
to generate an effective propagator equal to the identity) 
is an important benchmark for control |29| and will be 
explored in future work. 

The coherence times of C\ and C2 indicated in Fig. [S] 
(and the corresponding Lorentzian linewidths shown in 
Table Hljl demonstrate a ~ 50-fold increase in coherence 
time compared with the free-induction dephasing times 
(Table QJ. Further, the residual broadening is compara- 
ble to the strength of the C a -Cp dipolar interactions, as 
expected. The coherence time of the methylene carbon 
C m is clearly limited by the efficiency of proton decou- 
pling. We note that dephasing times > 21 ms were not 
explored due to the heating effects of high-power pro- 
ton decoupling on the room temperature probe/sample. 
Comparing the Hahn echo and dipolar/Hahn echo data, 
it is seen that Zeeman shift dispersion contributes ~ 90 
Hz (1.2 ppm) to the free-induction linewidths, whereas 
the intermolecular dipolar fields contribute ~ 25 Hz in 
the 1.6% dilution crystal. Sequences that efficiently av- 
erage all spi n interactions, such as the Cory 48-pulse 
sequence [26|, could be used to test the limits of line- 
narrowing. 
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TABLE II: Simulated and experimental dephasing rates for natural abundance spins in 1.6% isotopic dilution crystal, in terms 
of Lorentzian linewidths (all in Hz). Simulations I/II refer to the exclusion/inclusion of couplings between indistinguishable 
spins. 'Simulation IP values should be compared to the Hahn echo results, and 'Simulation P values should be compared to the 
dipolar+Hahn refocusing results. Under the latter sequence, effective linewidths are reduced by nearly two orders of magnitude 
compared to the free- induction decay spectrum, with the exception of C m , whose coherence time is limited by residual proton 
dephasing (see footnote a). 





Dipolar and static-field dephasing (1.6% crystal) 


Simulation I 


dipolar + 

Hahn 
refocusing 


Simulation II 


Hahn echo 


FID 


CI 


6.7 


~ 2 


20.1 


25 


133 


C2 


7.6 


~ 3 


25.1 


32 


122 


Cm 


13.9 


ll/20 a 


18.0 


25 


103 



a ll Hz and 20 Hz correspond to proton decoupling powers of 330 
kHz and 270 kHz, respectively. 



VI. DISCUSSION 

In this paper, we have established benchmark results 
and described a methodology for controlling solid-state 
NMR qubits upon which future work can build and im- 
prove. High-fidelity quantum control was demonstrated 
through the preparation of a pseudopure state. That re- 
sult, along with unitary simulations of the fidelity for the 
3-qubit controlled- (not(8not) gate, suggests that ensem- 
ble inhomogeneities involving linear Hamiltonian terms 
(Zeeman dispersion, RF amplitude) may be suppressed 
by strongly modulating RF pulses. Additionally, the 
method of strongly modulating pulses is well-suited to 
finding solutions for generating arbitrary quantum gates, 
approaching time-optimality, in systems with complex in- 
ternal Hamiltonians. Measurement of the triple-quantum 
dipolar dephasing rate was carried out, as an applica- 
tion of the pseudopure state protocol. Simulations of 
the intermolecular dipolar field spectrum predicted dipo- 
lar linewidths in reasonable agreement with experimental 
results obtained under appropriate refocusing sequences. 
Refocusing of the Zeeman and like-spin dipolar terms 
lead to a 50-fold increase in ensemble coherence times, 
and suggests that much narrower effective linewidths 
could be achieved using sequences that efficiently remove 
all spin interactions. 

These results demonstrate the potential for achieving 
universal quantum control in solid-state NMR systems 
with the high precision required for quantum computa- 
tion and other information processing tasks. Coupled 
with state purification techniques like dynamic nuclear 
polarization, such control should enable a powerful new 
QIP testbed reaching up to ~ 20 — 30 qubits @ . Fur- 
thermore, the solid-state setting is more general than the 
liquid-state, for example, by allowing one to couple the 
qubits to a bulk-spin heat bath as demonstrated by 



the selective polarization transfer that was employed to 
begin the pseudopure state protocol. Similar techniques 
will enable implementation of quantum error correction 
and heat-bath algorithmic cooling |lj] protocols, the lat- 
ter having been already demonstrated in the malonic acid 
system [lq . 

To further improve control, the strongly modulating 
pulse optimization procedure must be improved to yield 
simulated fidelities comparable to those achievable in 
the liquid state (that is, unitary fidelities in the range 
of 99 - 99.9% as opposed to - 95 - 99%). Moreover, 
the experimental system must be carefully engineered so 
that the implemented control fields are sufficiently faith- 
ful to the simulated pulses. Extending to larger num- 
bers of qubits (> 12) will likely require some form of hy- 
brid control that utilizes multiple-pulse, average Hamilto- 
nian techniques in combination with the numerical opti- 
mization methods discussed. Additionally, scalable con- 
trol methods that have been successfully demonstrated 
in the weak-coupling regime 0, || could potentially be 
merged with the strongly modulating pulse methodology 
to achieve a scalable pulsefinder for appropriately tai- 
lored systems that include strongly-coupled spins. 

In conclusion, we have characterized in detail a novel 
solid-state NMR system for quantum information pro- 
cessing, and used it to demonstrate high-fidelity quan- 
tum control. This lays a foundation for future studies in 
the malonic acid system, and in larger systems, desirably 
with high polarization. 
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